Nonperturbative Treatment of Double Compton Backscattering in Intense Laser Fields 
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The emission of a pair of entangled photons by an electron in an intense laser field can be described by 
two-photon transitions of laser-dressed, relativistic Dirac-Volkov states. In the limit of a small laser field inten- 
sity, the two-photon transition amplitude approaches the result predicted by double Compton scattering theory. 
Multi-exchange processes with the laser field, including a large number of exchanged laser photons, cannot be 
described without the fully relativistic Dirac-Volkov propagator. The nonperturbative treatment significantly 
alters theoretical predictions for future experiments of this kind. We quantify the degree of polarization cor- 
relation of the photons in the final state by employing the well-established concurrence as a measure of the 
entanglement. 
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Introduction. — In ordinary Compton scattering [1], a pho- 
ton is scattered inelastically by an electron. For photons with 
energy much less than the electron's rest mass, the quantum 
mechanical expression for the cross section agrees with the 
one obtained by classical electrodynamics. Nonlinear Comp- 
ton scattering is encountered when several photons from a 
strong laser beam are scattered by a free electron to produce a 
photon of different energy; this process has been calculated 
theoretically y, \M and successfully measured 14, bfl. Re- 
cently, there has been an increased interest in a different non- 
linear generalization of Compton scattering where a free elec- 
tron collides with a laser pulse and emits two photons at the 
same time. This process has no classical counterpart, and in- 
deed, as we will see, the two photons exhibit a paradigmatic 
quantum feature: namely, their polarizations are entangled. 
Properly optimized, two-photon emission from backscatter- 
ing of laser photons at an electron beam holds the promise of 
providing entangled light at much larger energy than conven- 
tionally used for quantum information purposes 

With relativistically strong lasers being available in many 
laboratories worldwide, the current record being a laser in- 
tensity of 10 22 W/cm 2 at the focus 01, the quest for observ- 
ing genuine laser-induced quantum effects in the relativis- 
tic regime continues. However, the peak field strengths are 
still orders of magnitudes below the quantum electrodynamic 
(QED) critical field E c = —m 2 /e = 10 16 V/cm for pair cre- 
ation (here m and e = —\e\ denote the mass and charge of 
the electron, respectively, and we use natural relativistic units 
c = h = 6q = 1). Two-photon emission by a laser-dressed 
electron via nonperturbative double Compton backscattering 
is a strong-field, relativistic quantum effect which could be 
observed without the additional complications connected with 
the ultra-relativistic particle beams necessary for laser-dressed 
pair creation fiHEfll. 

The theory of perturbative double Compton scattering, 
the reaction in which one photon scatters on an electron 
to produce two final photons was calculated by Mandl and 
Skyrme lUlll . recently reexamined in Jl2ll . and experimen- 
tally confirmed in Ed^]. The relevant Feynman diagrams 
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FIG. 1: The Feynman diagrams for perturbative (a), and nonpertur- 
bative (b) double Compton scattering. In (a), an electron with initial 
momentum pi absorbs one laser photon x, emits two photons with 
wave vectors kb and k c and ends with final momentum pf. Instead 
in (b), the laser-dressed initial state with average four-momentum 
qi, decays to the final state with average four-momentum qf under 
emission of two photons. In the intermediate state the momenta are 
labeled by pb and p c , respectively. The nonperturbative interaction 
with the laser field is depicted by dressing the electron lines with a 
zigzag line. 



are displayed in Fig. QIa). In Q1J 03], the simultaneous 
emission of two photons is interpreted in terms of the Un- 
ruh effect. Other two-photon processes that have been in- 
vestigated, both theoretically and experimentally are double 
bremsstrahlung tl8l [l9l |20|], two-photon synchrotron emis- 
sion [2jJ, |22|], and the total rate of two-photon emission in a 
crossed field l23ll . However, the generalization to nonpertur- 
bative double Compton scattering has not been recorded in the 
literature to the best of our knowledge. 

The purpose of this Letter is twofold: To show (i) that 
photon pairs with quantifiable entanglement can be produced 
from double Compton scattering in intense fields, and that ( ii) 
nonperturbative effects have to be incorporated to make reli- 
able predictions for a relativistically strong laser pulse. 

Nonperturbative QED formulation. — The interaction of an 
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electron with a laser field of arbitrary intensity can be treated 
in the formalism of strong- field QED (Furry picture), where 
the classical external field is included in the unperturbed 
Hamiltonian. For the present problem, the calculation of the 
amplitude of the process amounts to evaluating the Feynman 
diagrams shown in Fig. \Hb). The external electron lines 
are Volkov states exact solutions to the Dirac equation 
with an external laser field, (id — m — eA)^! = 0, and 



the propagator is the Dirac- Volkov propagator ||24j, 12511 . 
We denote four-vector scalar products of four-vectors q 
and p as q ■ p = q^p^ = q°p° — q ■ p, and the Dirac 
contraction is written as p = 7 • p. The laser four-vector 
potential A M = a M cos(x • x) propagates in the negative 
x 3 -direction with wave four- vector h and frequency u, and is 
linearly polarized in the x 1 -direction. We also introduce the 
intensity parameter £ = \e\^/\a 2 \/2/m, which can be used 
to classify the regime of relativistic laser-matter interaction: 
£ <C 1 corresponds to the perturbative regime, and £ > 1 
to the nonperturbative. To describe the emitted photons, 
we employ spherical coordinates so that the momenta read 

h,c = ^6,c(l, sin b , c cos Vv, sin 0b,c sin ^&,c cos 0b,c), 
with being the frequency. To study the po- 

larization correlation, we also need a basis for the 

we use 
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polarization vectors of the photon pair: 

cos Vv,cos #6, c sin i/j b ,< 
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generic polarization vector of the photons is 

66,c = (cieJ jC + c 2 e 2 hc )/^\c 1 \ 2 + |c 2 
constants c\^- The initial electron is assumed to propagate 



given by 
for some complex 



in the x 6 -direction with four-momentum pi 



colliding head-on with the laser pulse. The average, or quasi 
momentum J26ll of the electron immersed in the laser wave is 
given by q { = pi - xe 2 a 2 /(4x • p^ = (Q*, qi), with average 
mass m* = \/gf, and the corresponding quantities for the 
final electron are labeled by pf, qf. 

Having fixed the notation, we proceed to calculate the scat- 
tering amplitude S. The calculation follows the usual steps of 
laser-dressed QED |25|], and we present only the final result, 



S = i 
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The energy-momentum conserving delta function contains 
the integer n, which is the net number of photons absorbed 
during the entire collision. The second index of summation s, 
which appears in the propagator momenta p^ c , is the net num- 
ber of photons exchanged before emitting the second photon 
[k b or k c , depending on the diagram, see Fig.[TJb)]. The am- 
plitude 0} is gauge invariant under e b , c — > e b ^ c + Ak b , c , where 
A is an arbitrary constant. Another important aspect of the 
amplitude S is the possibility for the propagator momenta to 
reach the laser-dressed mass shell p 2 ch = m 2 , which indicates 
the split up of the process into two sequential single Compton 
scattering events |23|]. At such a resonance, where the ma- 
trix element formally diverges, S may be rendered finite by 
including a small, imaginary correction to the laser-dressed 
electron mass and energy HHE], or alternatively be regular- 
ized with an external parameter such as the laser pulse length 
or a finite detector resolution. In the following, we will always 
consider parameter regions such that the sequential Compton 
scattering cascade is forbidden by energy-momentum conser- 
vation or is exponentially suppressed by a large-order Bessel 
function. This selection is in accordance with planned exper- 
iments recently discussed in Refs. HO, EH]. In order to facil- 
itate the detection of the rather weak two-photon signal, the 
measurement should be done in energy and angular regions 
where the the single Compton scattering process is strongly 
suppressed. 

In the following we evaluate the differential rate 
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where T is the long observation time. Integrating over the 
final electron momentum qf and the final photon energy uu c , 
we end up with the rate dW / duobdVtbdVt c , differential in the 
directions dftb,c = dcos 6b, c dipb,c of the two photons and in 
dujb- The photon energy u c is given by 
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Here V is the quantization volume, M N 
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-r^- 2 — k, and A{ AT is a generalized Bessel function 

(27|], A?* N = J 27r cos fc 6> c iAr6>-ifa,-a^sin6>+if/3,-/3^sin26>^ 

, h G {0, 1, 2}, with the arguments aj = ea ■ Pj/(x ■ pj), 
(3j = e 2 a 2 /(8x • pj), j G {z,/, 6, c}. The spinors uij are 
normalized according to u\ fj°Ui j = 1. 



where k c /uo c is independent of uo c . The last line in Eq. ® 
holds if <C -gr w Q b w Q c <c 1, which is the parame- 
ter regime on which we will concentrate (backscattering ge- 
ometry with relativistic electron energy). Moreover, Eq. ([3]) 
implies that the sum of the two photon energies is limited by 
uj h + u c < 4nco(Ei/m) 2 /(l + £ 2 ). 

Calculated differential rate. — In Fig. [2 we show the differ- 
ential rate in the laboratory frame, for a specific set of parame- 
ters, and compare to the corresponding rate obtained from the 
perturbative formula llll LL2] which includes only one inter- 
action with the laser field [see Fig.QJa)]. We have checked 
that the expression © agrees with the one obtained in Willi in 
the limit of small laser intensities. Since we take the initial 
electron to be relativistic, Ei = 10 3 m, it will emit mainly 
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FIG. 2: (color). Comparison of the perturbative and nonperturbative 
approach. Shown is the laboratory frame rate dW / dubdQbdQ c , in 
units of s _1 sr _2 MeV _1 , for the nonperturbative [(a), (c)] and per- 
turbative [(b), (d)] case. The parameters employed are Ei = 10 3 m, 
ou = 2.5 eV, £ = l,u b = l MeV, and b = 6 C = 10" 3 . The photon 
polarizations are given by 6b, c = e\,c m ( a ) an ^ 3 ), and 6b, c = ^l,c 
in (c) and (d). 



in the forward direction, 9b, c ~ m/Ei. The laser parameters 
= 2.5 eV and £ = 1 corresponds to an optical laser with 
intensity 5.5 x 10 18 W/cm 2 . Since the quantum parameter 
X — iVi ' x/m 2 1 2611 is small (= 10 -2 ) here, spin effects 
are marginal and we therefore average (sum) over the initial 
(final) spin of the electron. The small value of x a l so P er_ 
mits us to neglect effects arising from electron-positron pair 
creation, since the e + -e~ production rates are exponentially 
suppressed. For the parameters used in Fig. O up to n = 20 
laser photons participate, so that u c < 60 MeV according to 
Eq. ([3]). The results from Fig. [2] suggest that the differential 
rate varies strongly as a function of the angles and polariza- 
tion. It becomes clear that to interpret data from planned ex- 
periments of this kind 13011 , the nonperturbative formula © 
has to be used. 

To answer the question whether the integrated nonpertur- 
bative rate differ significantly from the one predicted by the 
usual double Compton scattering formula, we show in Fig. [3] 
the differential rate dW/dO c , integrated over the azimuth an- 
gles ip^ c , the polar angle Ob of one of the final photons and 
the energy cob, and summed over the final photon polariza- 
tions. The energy integration was limited to the interval be- 
tween 1 keV and 1 MeV to avoid the infrared divergence at 
ujb — > l32ll and cascade contributions for larger and 
for the same reason the integration over Ob was performed 
over the interval (0, 1.5 x 10 -3 ) radians. Restricting the final 
phase space in this way ensures that contributions from sin- 
gle Compton scattering are negligible; at polar angles smaller 
than 1.5 x 10 -3 radians all harmonics occur at energies larger 
than 1 MeV Integrating the nonperturbative curve in Fig. [3j 
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FIG. 3: (color online). The integrated and polarization- summed lab- 
frame rate dW/d0 c = sin c dW /d cos C , differential only in the 
angle 6 C . Shown are results for the perturbative and the nonperturba- 
tive calculations, for Ei — 10 3 m, 00 = 2.5 eV, £ = 1. 



we obtain a total rate in the laboratory frame of W = 3.5 x 10 T 
s _1 . For the perturbative curve, we get Wp er t = 2.5 x 10 7 s _1 , 
from which we gather that even for the integrated rate, the 
nonperturbative corrections are significant. The obtained two- 
photon rate should be compared to the total rate of nonlinear 
single Compton scattering l26h . which amounts to 3 x 10 13 
s _1 for the same parameters as in Fig.[3j Employing an elec- 
tron beam with 10 9 electrons per bunch, a laser pulse of dura- 
tion 100 fs, photon energy uj = 2.5 eV, intensity 5.5 x 10 18 
W/cm 2 (corresponding to £ = 1) 113 311 . and assuming perfect 
transverse overlap of the two pulses, we estimate that about 
2 x 10 3 photon pairs per shot may be expected. 

Entanglement. — Having investigated the differential and to- 
tal photon pair production rate, we now turn to the interesting 
question of the quantum mechanical correlation between the 
final state photons. To quantify the degree of polarization 
entanglement, we employ the well-established concurrence 
C(pf) D4f] as an entanglement measure. Assuming an un- 
polarized initial electron and unobserved final spin, we trace 
out the spin polarizations of the initial and final electron and 
calculate the 4 x 4 final density matrix p f of the polarizations 
\b, c £ {1,2} of the two emitted photons. Then, C(pf) is 
given by 

( 4 ) 



CO*,) = max(0, VCi - VC2 - VCs 

where the Q 's are the eigenvalues, in descending order, of the 
matrix Q = pf(cr y <8> a y )p^(a y 0a y ), where a y is the second 
Pauli matrix. For a maximally entangled state, C(pf) = 1, 
and for a non-entangled state C(pf) = 0. We note that the 
concurrence has recently been used to study correlation in the 
two-photon decay of a bound state 1135m . For the present case, 
the final density matrix p , can be computed from the normal- 
ized matrix element (Q]). Writing S = S rurf (A5, A c ), where Ti 
(r,) denotes the spin polarization of the initial (final) electron, 
we have for the matrix elements of pf, 

(\b, AcIp/IAj, a' c ) = — s rii r f (\b, K)s*i,r f (Ki K)- 

ri,r f 

(5) 

Here N is a normalization constant, which can be found by re- 
quiring Tr pf = 1. The concurrence, as defined in Eq. ©, is 
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FIG. 4: (color). In panel (a) we show the concurrence C(pf) [see 
Eq. ©], using the nonperturbative expression for the matrix ele- 
ment, which should be compared to panel (c), where the concurrence 
is shown for the perturbative case. The parameters employed are 
Ei 10 3 ra, oo 2.5 eV, f = 1, uj h = 1 MeV, and ip b = 2p c = 0. 
As a reference, we also display in the right column the logarithm of 
the polarization- summed differential rate \og 10 (dW /dcubdQbdQ c ), 
in units of s _1 sr _2 MeV _1 , for the nonperturbative [(b)] and pertur- 
bative [(d)] case. 

a gauge invariant quantity, furthermore it does not depend on 
the basis used to describe the polarization of the photons k^^. 
C(pf) depends sensitively on the energy and the directions 
of the emitted photons. One example of the fully differential 
concurrence is displayed in Fig. [U which shows the neces- 
sity of the nonperturbative formalism to predict the degree of 
entanglement. 

Conclusions. — We have studied the process of nonperturba- 
tive two-photon decay of a laser-dressed electron. Our results 
significantly alter the theoretical predictions as compared to 
a perturbative treatment of the laser; they lead to novel fea- 
tures in the angular and integral characteristics, which could 
be resolved using presently available intense laser facilities. 
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